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Abstract 

D , 
,i^ . We develop an algorithm which can be used to exclude the existence of classical 

breathers (periodic finite energy solutions) in scalar field theories, and apply it to 
S^ . several cases of interest. In particular, the technique is used to show that a pair of 

j^ . potentially periodic solutions of the 3+1 Sine-Gordon Lagrangian, found numerically 

in earlier work H, are not breathers. These "pseudo-breather states" do have a sig- 
nature in our method, which we suggest can be used to find similar quasi-bound state 
configurations in other theories. We also discuss the results of our algorithm when 
applied to the 1+1 Sine-Gordon model (which exhibits a well-known set of breathers), 
and cp'^ theory. 
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1 Introduction 

The leading order chiral lagrangian description of pions, when evaluated on an isospin- 
polarized field configuration, is identical to the 3+1 Sine-Gordon lagrangian 

C = ]^d^<Pd^(t> + Flm'cos{(t>lF.) ■ (1) 

The motivation for considering isospin-polarization of pion fields arises from the possibility 
of producing disoriented chiral condensate (DCC) |^ configurations in heavy ion collisions. 
In previous work [Q, we discussed the numerical properties of a pair of spherically symmetric, 
long lived classical solutions which arise from this lagrangian (related studies appear in Q 
and 0]). The evolution of the two configurations is nearly periodic, with energy localized 
for hundreds of oscillation cycles. The properties of the lower energy solution are such that 
both the quantum corrections and higher derivative terms of the Chiral Lagrangian should 
be negligible. Hence it is likely to correspond to a real physical pion configuration in QCD 
(pionic breather state, or PBS), which could in principle be produced in heavy ion collisions. 
Although many of its physical properties have been studied in , some theoretical questions 
remained unanswered, such as whether truly periodic, finite energy solution might exist in 
the 3+1 Sine-Gordon model. 

In this paper we show that the two PBS solutions do not correspond to true breather 
solutions, where a breather is defined to be a periodic, finite energy classical solution. Our 
method can be summarized as follows. First, periodic solutions (f){r,t) are classified in the 
large r region based on their asymptotic behavior, in terms of a finite number of parameters. 
Then a numerical algorithm is used to evolve the configuration using the equations of motion 
from large r to r — > 0. Most choices of asymptotic behavior lead to singular behavior at 
r = 0, and hence do not correspond to finite energy solutions. Using the criteria of finite 
energy, we can search the parameter space in a systematic manner. In section 2, we review 
the numerical properties of the two PBS solutions. In section 3, we give an overview of 
our algorithm. In section 3.1, we analyze the free Klein-Gordon equation and its full set of 
periodic solutions. In section 3.2, we use the solutions of section 3.1 as asymptotic forms for 
the full Sine-Gordon theory to search for breather solutions, and present the results. Section 
4 contains a proof that interactions do not enlarge the set of asymptotic solutions beyond 
those of the free Klein-Gordon equation. That is, there is a one-to-one mapping between 
asymptotic behaviors of the free and interacting equations of motion. In section 5 we analyze 
the 1+1 Sine-Gordon and (p^ theories, finding that the former theory contains exact breather 
solutions while the latter does not. In Section 6 we offer concluding remarks. 



2 A review of PBS properties 
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Figure 1: Evolution of the lower energy PBS 



Both PBSs are spherically symmetric. Imposing spherical symmetry, the equation of 

motion is 

0_0"_^0' = _sin0 , (2) 

r 

where / = f and /' = £. 

Figures |I] and ^ show the numerically generated evolution of the two breather candidates, 
exhibiting the 0(r, t) vs. r dependence at different times. The nearly periodic evolution 
continues virtually unchanged for hundreds of cycles but eventually breaks apart, dissipating 
energy to radial infinity. 

The frequency and approximate lifetime of the lower energy PBS are 0.1 cycles/ (fm/c) 
and 850fm/c. For the higher energy PBS, they are 0.0619 cycles/(fm/c) and 4000fm/c. 
Additional searches have not indicated a third PBS state. 

By an appropriate scaling of and x'^, the Sine-Gordon lagrangian (equation |l]) can be 
written as 

(3) 



C = - d^cj) df^(j) + cos . 

In these units, the numerically calculated frequencies of the PBSs are uj = 0.9 ± 0.05 and 
cj = 0.58 ± 0.02. 
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Figure 2: Evolution of the higher energy PBS 

3 The Algorithm 

For finite energy solutions, ^^Oasr^ooQ. As 0^0, nonlinearities in the equation of 
motion become negligible, hence the r -^ oo behavior of any breather solution must approach 
a periodic solution to the noninteracting Klein-Gordon theory. 

The strategy of the following algorithm is as follows. First, we solve for the most general 
periodic, finite energy solutions of the Klein-Gordon equation. These solutions are used as 
large r asymptotic forms and are evolved inward to the origin using the full Sine-Gordon 
equation of motion. Only solutions which are finite and differentiable at the origin can 
correspond to true breather states. Since the space of asymptotic forms is parameterized by 
a finite number of variables, a comprehensive search can be initiated which in the end can 
rule out the existence of breathers within a certain frequency range. 



^Sine-Gordon finite energy solutions could converge to any multiple of 27r, but these solutions are equiv- 
alent by a translation in 0-space to those with — > 0. 



3.1 Asymptotic Forms 

It can be shown that the most general periodic, spherically symmetric classical solutions to 
the free Klein-Gordon equation are linear combinations of 

(pUr, t) = - sin(fcr - 9k) sm{ujt - 9^) , (4) 

r 

where 

w'^ = k'^ + m^ . (5) 

If TTi = 0, these solutions have infinite energy for all cj, and no breathers exist. For nonzero 

m, the set of solutions for {) < uj < m have imaginary k, corresponding to the finite energy 

solutions of the form [| 

4>uiir,t) = - exp(—kr) siniut — 9) , (6) 

r 

where 

(Since these are the only types of solutions we are interested in, we will drop the tilde from 
now on). The most general finite energy periodic solution is therefore 

= - ^ ajexp{—kjr)sm{juft — 9j) , (8) 



where uif is the fundamental frequency in the Fourier series, and kj = Jm? — {oojy. Note 
that in the case that 0.5m < oOf < m, this sum has only one term 

= - aiexp(— fcir) sin(ci;/t) , (9) 

hence the full set of solutions is parameterized by two variables, ai and Uf. The frequencies 
of the two numerical PBS solutions fall well within this frequency range. 

In the remainder of this section we will work in units with m = 1. 



3.2 The Search for Breathers 

Using the asymptotic form in equation ^ as r ^ oo boundary conditions, the Sine-Gordon 
equations of motion can be integrated inward to r = (see appendix A for computational 
details). In general, solutions obtained in this way are neither finite nor differentiable at 

^These are also not breathers, because (/)[^(0,i) ^ 0, and hence there are no true breathers of the free 
Klein-Gordon equation. 



the origin, and therefore do not correspond to physical breather states. For small r, the 
numerically obtained solutions can be fit to the form 

^ + B{t) , (10) 

r 

and the integral 

f27r/ujf 

S= dt \A{t)\ (11) 

Jo 

can be used as a measure of how singular each solution is at the origin^. A similar measure 

of differentiability could be defined (for instance by fitting to a form A(t) r + b(t)), but as 

will soon be shown, in the present case not even the singularity condition is fulfilled, so it is 

not necessary. 

The size of the singularity S is calculated for various values of ai and u. To fill ai — u 
space as uniformly as possible values are randomly chosen in the range 0.5 < Uf < 1 and 
Oi > 0. In practice an upper limit on ai has to be set, so a value is chosen which is much 
larger than that of the two known pseudo-breather states. Also, the lower limit of ai had 
to be raised slightly (from to 1.2), because ai —>■ corresponds to the trivial solution 
0(r, t) = 0, which lowers the value of S to zero across all values of uj. 

Figure ^ shows a plot of S vs. ujf (the value of ai is not shown, so the figure can be 
thought of as a transparent projection plot). As can be seen, this plot has three prominent 
minima. Two occur at the frequencies of the previously discovered PBSs, u ~ 0.560 and 
u ~ 0.865. Neither minima goes all the way to zero, and hence they do not represent exact 
breathers. The third minima at a; = 1 corresponds to the trivial solution 

0(r,t) = Asm{mt-e) , (12) 

with infinitesimal A. 

As further evidence supporting our hypothesis that these minima correspond to the 
previously discovered PBS states, figure | shows the shape of a 0(r, t) calculated near the 
minima at u ^ 0.560. Notice the similarity with figure ^ with the exception of r ^ where 
figure ^ becomes singular. 

We now define the solutions corresponding to the minima of the S vs. Uf curve as 
'pseudo-breather' states (a useful reinterpretation of the acronym 'PES', which previously 



■^Although S is obtained by fitting to a 1/r form, a nonzero value of S can be shown to also indicate 
singularities of the type </> ~ 1/r" with n > 1. A first derivative fit near the origin would be even more 
powerful (capable of fitting derivative and logarithmic singularities), but tends to yield ever increasing values 
of S as the fit is done closer to the origin, which is less manageable. It will therefore be avoided until cases 
where it is necessary, such as the 1+1 Sine-Gordon equation discussed later in this paper. 
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Figure 3: The singularity S vs. oscillation frequency uj 

referred to 'Pionic Breather State'). Although not true breathers, pseudo-breathers can often 
exhibit very interesting behavior. As was discussed in 0, the two PBS configurations arise 
readily from multi-domain DCC configurations, have characteristic energies and frequencies, 
and might extend the lifetime of classical pion configurations by more than an order of 
magnitude. 



4 A More Formal Justification 



At this point questions still remain about the comprehensive nature of the solutions obtained 
in the previous section. It is unclear how the addition of interactions might enlarge the full 
set of periodic solutions. Perhaps each asymptotic form corresponds to multiple solutions, 
and integrating inward only uncovers one. In this section we reanalyze the equations of 
motion by expanding in a power series. Without solving the resulting complicated set of 
equations, we study its general nature to answer the question just presented, and therefore 
verify that the previously described algorithm is valid. 
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Figure 4: Evolution of 0(r, t) at the u ~ 0.560 minima of figure ^ 

In section 4.1, we work out a general form of the solutions. In section 4.2, we use this 
general form to solve the free Klein-Gordon equations as a test case for the full theory, which 
we finish with in section 4.3. 



4.1 A Series Form for the Solutions 

First we show that the most general periodic solution with frequency 0.5 < tu < 1 is of the 
form Fl 

p—kr oo oo 

0(r,t) = a sm{ujt) — — + Y. J2 ^ii.iaW 



-J2kr 



ji=2J2=2 

where the ^jual^) ^^^ unspecified functions, and k'^ + uP' 
Parameterizing the interaction strength as A, 

£ = \d^^d^,^ - \m^^ + A r (0) 
we can now write the solution as a series in A, 

oo 

0(r,t) = ^A>„(r,t) . 

n.=0 



rJi 



2 



(13) 



(14) 



(15) 



"^In many particular cases such as the Sine-Gordon or (/)^ lagrangians, the sum in this general form starts 
at ji = 3 and 22 = 3. This is however a particular case of the form as given, and therefore need not be 
considered separately. 



We will inductively show that each term is of the form of equation ^ 



Substituting |T5| into the equations of motion, yields 

[9' + m2] 0„(t,r) = P„[0o(r,t), 0i(r, t), .... 0„_i(r, t)] 



(16) 



where Vn are polynomials expressions related to the Klein-Gordon interaction term in a 
complicated way. The equation for n = is just the Klein-Gordon equation of motion 







hence its solution is the usual one given by equation P, 



;i7) 



a sm{ujt) 



-kr 



This is also the first term in equation 13 



Iteratively substituting in the expressions for 00; 0i; ••••0n-i, we can solve for 0„. 

To show that higher 0„'s are in the form of the remaining term in equation |l^, first 
note that since V is polynomial whose lowest order term is at least quadratic, inserting 
00, 01, ••••0n-i5 yields something of the form 

-J2kr 



oo oo 



>^)=E H^n.M 



ji=2J2=2 



fn 



(19) 



(this differs from equation 13 only in that it is missing the first term). Equation 16 



IS a 



second order linear inhomogeneous differential equation whose most general solution is 

g-fcr 

0„(r,t) = a sin(wt) + Q(r,t) , (20) 

r 

where a is an undetermined variable and Q(r, t) is any particular solution. It can be shown 
that [d"^ + rn?] Q spans the space given by equation [1^ if Q is of the form in equation |T3|. 
A Q in the desired form can therefore be chosen to satisfy the differential equation, and by 



equation |2y, 0„ is also in this form. 



4.2 The Klein-Gordon Case 



Before we work with the full equation of motion, we will use the form of equation O in the 
free Klein-Gordon Lagrangian to derive the full solution for a given large r behavior of (in 
the form of equation P), and show that this is unique. Although this result is trivial, it sets 
up a formalism which can be naturally extended for the next section when we will repeat 
this for the full theory including interactions. 
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Substituting the form of equation ^ into the Klein-Gordon equation of motion yields 
the decoupled set of equations 

^,.2(i) + Hi k' + m']A2At) = , (21) 

Asj, (t) + HI k' + m']A3,,, (t) = Aj2kA2,,, (t) , (22) 

and then for ji > 3, 

An,hit) + Hi k^ + "^''\K,n{^) = 2j2fc(ji - 1)A,,^,,jM + JiUi " l)Ai-2,,,(t) • (23) 
These are each second order linear inhomogeneous differential equations with general solution 

Aji,j2 (t) = (^h,h sin(wi2^ - Qn,n) + ^ii ,i2 (*) , (24) 

where a^^^^ and Qj^^j^ are undetermined constants, Pj^^j^it) is any particular solution, and 

uj],^{k22f = m' . (25) 

Note however that each a^^^j^ is forced to be zero. The ^vniujj^t — ^jija) is an on-resonance 
driving term for the next differential equation (with ji incremented by one), and the asso- 
ciated solution is not periodic. All freedom in the choice of a^^jj is lost. The driving terms 
for each differential equation iteratively become zero and therefore each Pj^^j^ii) = 0. 

The most general solution with the given lowest order asymptotic term is therefore 

4> = - a exp(— /cr) sin(u;t) . (26) 

r 

4.3 Adding Interactions 

Now we repeat the analysis in the last section for the case where there is an interaction 
term present, and arrive at our final conclusion for this section, that the solution to the full 
equation (interaction present) is completely determined by the frequency and size of the first 
order asymptotic piece that is a solution to the free Klein-Gordon theory. 



Repeating the same analysis with an interaction term changes equation |21| - p3| to 

^2,n{t) + I-jI e + m^]A2,j,{t)=Vij,[A,,uA2,2,.-] , (27) 

AsAt) + Hi k" + nn?]A^^j,{t) = 4j2kA2j,{t) + V2,j,[Ai,i, ^2,2, ...•] , (28) 

and then for ji > 3, 

^n,hit) + Hi k' + ^Mjuhit) = 2j2A:(ji - l)A,-,_i,,,(t) + ji(ji - l)A,,_2,,,(t) + (29) 

^iij2[^l,b^2,2,----] • 

10 



The actual form of Vj^j^ is related to the Klein-Gordon interaction term in a complicated 
way, but two things are clear. First, since the interaction contains no linear terms, Vj^j^ 
is at least quadratic in the Aj^jj's, and can not explicitly cancel the linear terms already 



present from equations ^l| - |23| . Second, each "Pjija o^^Y contains terms A--^-j^{t) with ji < ji 
and J2 < J2- This allows us to solve iteratively beginning with the lower order solutions, 
yielding once again simple second order linear inhomogeneous differential equations. 

Unlike in the noninteracting case, after inserting in the lower order solutions Aj-^^,j,^{t), 
Vj-^j^ may have a term in it proportional to the on-resonance driving piece, sm{uj^t — Oj-^j^). 
This doesn't add any freedom in our choices of aj^,J2 however, in that now these values must 
be chosen to cancel the resonant pieces. This constraint completely specifies the form of 
each Aj-^j2{t), hence there is no additional freedom in the choice of solutions. 

The final solution is therefore unique given its lowest order large r asymptotic behavior. 



5 Other Examples 
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Figure 5: S vs. uj for the 1+1 Sine-Gordon Lagrangian 

We end by applying our method on two other examples, the 1+1 Sine-Gordon equa- 
tion, and 0^ theory. The 1+1 Sine-Gordon Lagrangian is a well studied model with known 
breathers, and is thus used as a test of our algorithm. 

Minor modifications had to be made to work in 1 + 1 dimensions. The asymptotic form 
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Figure 6: S vs. oj for the 3+1 0'^ theory 



in equation pi was changed to 



Oi exp(— /cir) smiujft) . 



(30) 



Also, in 1+1 dimensions, (f){r,t) does not readily become singular at r = 0, hence we tested 
the differentiability of at zero. 

The S vs. u plots for each Lagrangian appear in figures ^ and p. In the 1+1 Sine- 
Gordon plot the values of S go to zero across all values of u, indicating a class of breathers 
as expected. The 0"^ plot stays positive for all values of u, hence in 3+1 dimensional 0^ 
theory, there are no breathers in the frequency range 0.5 < a; < 1. 



6 Conclusions 

In this paper we have developed a numerical search algorithm for spherically symmetric 
breathers in spin zero theories. This algorithm was used to show that breather-like states 
found in a previous work were not actual breathers, but rather a new type of object, a 
'pseudo-breather', which, while long-lived, is ultimately unstable. We showed that the 
asymptotic behavior of periodic solutions which might have finite energy can be parame- 
terized in terms of a finite number of parameters. In the case of interest, 0.5m < ujf < m, 
there are only two parameters: frequency uj and magnitude ai. An extensive search within 
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the Sine-Gordon model failed to uncover any values of cu/, ai for which the solution is non- 
singular at the origin. The results for this model are shown in figure ^ 

We have, as of yet, left the region < Uf < 0.5m unexplored {uj > m is automatically 
excluded by large r energy considerations). It is straightforward in principal to extend the 
algorithm to arbitrarily small u. This is done by generalizing the form of the asymptotic 
solutions. For instance, in the region 0.25m < uj < 0.5m, the asymptotic form becomes 

= - Oi exp(/i;ir) sin(u;/t) H — a2exp(A;2^) sin(2c<jjt — 6*2) . (31) 

There is a numerical trade-off, in that to cover smaller regions of uj the number of parameters 
to describe the asymptotic solution increases, eventually to infinity as cj/ ^ 0. 

The condition of spherical symmetry can be relaxed slightly by adding higher partial 
wave terms 

^{x,t) = 0°(r,t) Y,^{9,^) + <P^\r,t) Ff ^(^,0) + 0?(r,t) Y^{e,<P) + <p\{r,t) r/(^,0) + .... 

(32) 
Inserting this into the equation of motion yields a coupled set of equations 

c - c" - - 0r + ^^^0r = -v I E ^rc I • (33) 

^ ^ \l,m j 

Again, the number of variables needed to parameterize the space of asymptotic forms in- 
creases, thus demanding more computational power for each spherical harmonic added. 

Although not discussed here, our algorithm can be modified to work with higher spin 
theories. The condition of spherical symmetry greatly restricts the asymptotic forms, so in 
many cases not much more computational power is needed. 

Two shortcomings of the analysis in this paper are as follows. First, we can not guarantee 
that our upper cutoff in a\ hasn't hidden a true breather, and second, it is possible (as often 
with a numerical analysis) that some exceptional behavior is concealed in the gaps between 
calculated values in the S vs. uj^ graph. Due to the large range of points and small spacing in 
figure 0, this seems unlikely. However, it prevents the result from being completely rigorous. 

The authors would like to thank R. Seals, C Buragohain, Alan Chodos, M. Gokgedag, 
Michael Ibrahim, and Vincent Moncrief for useful discussions and comments. This work was 
supported in part under DOE contracts DE-AC02-ERU3075 and DE-FG06-85ER40224. 
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7 Appendix A : Details of Numerical Evolution 

The spherically symmetric 3+1 Sine-Gordon equation of motion is 

0(r, t) - 0"(r, t) - -0'(r, t) = - sin(0) . (34) 

r 



Discretizing, and using the approximations 

dx 2Aa; 

and 

d^f ^^ /(x,+i) - 2/(x,) + /(a:,_i) 

dx'^ Ax^ 

the equations of motion can be written 

1 



(35) 
(36) 



[ri-i,tj) 



i- 1 



2<P{ri,t,)i-<P{ri+^,tj){i + l) + /\x' I f , (37) 



where 

/ = ^in^t^^^) + ^ir.^^^~^)-^<Pir,^h) + ^i^(^(,^^ ,^.)) _ (3g) 

Spatial boundary conditions are fit to an asymptotic form at r — >■ oo. In practice this 
boundary fit is done at a large but finite value r = tq. Boundary conditions in the t direction 
are of course periodic. The values ro, Ax and At were varied until the resulting evolution 
(and hence the shape of figure |^) converged. 
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